Exchange assisted tunneling and positron annihilation on inner atomic shells 
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It is known for a long time that the long range asymptotic behavior of the Hartree-Fock orbitals 
is different from that of the orbitals in the local potential. However, there is no consensus about 
observable physical effects associated with this asymptotics. Here we argue that weaker decrease of 
the Hartree-Fock orbitals at large distances is responsible for the positron annihilation on the inner 
shell electrons, which is observed experimentally. 
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Long distance behaviour of electron orbitals 

It is sometimes assumed that exchange interaction be- 
tween electrons is important only at short distances and 
can be neglected when one of the electrons is far from the 
origin. It is easy to see that this is not true if we consider 
long range asymptotics of an inner shell electron [IH3] • At 
large distances the exchange term for all occupied orbitals 
has the form r~ u (j) Vl where <p v is the outermost orbital 
with the highest energy e v . The power v depends on the 
leading multipolarity of the exchange interaction. The 
monopole component does not contribute to the asymp- 
totics due to the orthogonality of the orbitals and for the 
dipole component of the exchange interaction v = 2. 

This fact was first realized by Handy et al. PQ, who 
showed that generally speaking asymptotic behavior of 
all Hartree-Fock orbitals is given by the exponential with 
the highest energy e v (or the smallest binding energy): 

4>f \r»r v ~ r-"« exp (-V=2i» , (1) 

where v% is specific for each orbital and r v is the radius 
of the outermost atomic orbital. We use atomic units 
h = m e = \e\ = 1 unless stated otherwise. Because the 
monopole component of the exchange interaction does 
not contribute to the asymptotic behaviour, the expres- 
sion ([I]) does not apply to the systems where only s- 
orbitals are occupied. 

Later the asymptotic behaviour of the electron orbitals 
was reanalyzed by many authors [2HB]. Dzuba et al. [2] 
showed that the asymptote ([I]) holds for the relativistic 
Hatree-Fock-Dirac equations as well. The role of cor- 
relations were studied by Morrell et al. [J| and Katriel 
and Davidson [5], who demonstrated that Eq. holds 
also for the natural orbitals <pf at with nonzero occupation 
numbers. Natural orbitals are the eigenfunctions of the 
one electron density matrix p(x',x), where x = r, a and 
they can be found for any many electron wavefunction. 

Natural orbitals are mostly used in the context of 
the configuration interaction approach. For the core- 
valence correlations the many-body perturbation theory 



(MBPT) is usually more efficient. Within MBPT ap- 
proach Flambaum [3] showed that Eq. |l]) also holds for 
the Brueckber orbitals. The latter are solutions of the one 
particle equation with the correlation potential E added 
to the Hartree-Fock potential. The nonlocal part of the 
potential E decreases faster than the exchange potential 
and therefore does not change the long range behaviour 
of the Brueckner orbitals. 

In the solid state physics the long range exchange 
induced interaction is well known. For example, 
the Ruderman-Kittcl-Kasuya-Yasida (RKKY) exchange- 
induced spin-spin interaction is responsible for the mag- 
netic ordering in metal alloys [7HS]- The long range in- 
teraction in one- two- and three-dimensional systems was 
also considered in [3]. 

In spite of the results cited above and numerous other 
studies it is sometimes assumed that asymptotic behavior 
is an artifact of the used approximations and physi- 
cally observable one particle asymptote should depend on 
the energy Si of the orbital in question (see, for example, 
Ref. [TD]): 

<t>i\r->r v ~ r~ v * exp (-\f-2elr) . (2) 

Such asymptotic behavior follows, for example, from the 
analytical continuation to the negative energies of the 
wave function in the scattering theory [5]: 
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^ N _ _f0£ ei (fc iri+ ^iogn) x ^^ r - 1 (2...AT) 5 (3) 

where So,n is the matrix element of the S matrix and 
^f 1 is the wavefunction of the ion. However it is easy to 
present an example where Eq. ([2| is incorrect: the double 
well potential with one electron below and another above 
the barrier (see Appendix |A| and also Ref. [6]). 

To get asymptote ^ it is sufficient to assume that 
when the first electron is far from the origin, r± 3> 1, the 
N particle wavefunction ^ N can be factorized as 

9 N \ ri>1 =<P i (l)x*?- l (2...N). (4) 
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The TV-particle Hamiltonian can be written in a form: 

h(l) + h{2) ■ ■ ■ + h(N) + — 

k<i ™ 
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(5) 



In the last line we left two first terms in the expan- 
sion of 1/ru and skipped for simplicity the angular fac- 
tors P^cos^u). Substituting Eqs. Q and ([5| in the 
Schrodinger equation H Ni & N = E N "^ and integrating 
over the coordinates 2. . . N with the function ty* Ar_1 we 
get: 
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Due to the parity selection rule the second term of the 
multipolar expansion from ^ turns to zero here. The 
solution of this equation satisfies Eq. ^ with the energy 

e t =E N - if" 1 . (7) 

The obvious problem with the anzatz Q is that it is 
not antisymmetric in permutations. It also does not ac- 
count for the correlations between the outgoing electron 
and the electrons of the remaining ion. Let us see what 
happens if we substitute Q with the wave function of the 
general form *f? N . Following [5] we can project iV-particle 
wave function on the eigenstates of the (N — l)-particle 
ion: 

|* Ar (l...iV)>=^/ l (l)|*f- 1 (2...7V)), (8a) 
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Ml) = (^- 1 (2 . . . N)\^ N (1 . . . N)) 2 ... N . (8b) 



Expansion (8a) is valid at all distances and functions fi 



play the role of one-particle orbitals. They do not form 
an orthonormal basis set. In general these functions are: 
(i) not orthogonal to each other; (ii) normalized so that 
12i(fi\fi) — 1. (iii) n °t described by a definite angular 
momentum j.;. It follows from Eq. ( 8a ) that each function 
fi has several angular components: 



J+Ji 

fi(r) = ki,m(r) 



(9) 



where m — M — Mi and J, M and Ji , Mi are the angular- 
quantum numbers of the initial atom and the final ion 
respectively. 

It is easy to see that the amplitudes fi satisfy following 
equations [5]: 
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(10b) 



where p i fe ~ (r',r) is the transition density matrix of the 
ion. 

If we look for the solution of this system at large dis- 
tances it is useful to single out from the functions W^i 
the term with zero multipolarity and add it to the left 
hand side: 




fi(ri) 



r\ *■ — ' 



(11) 



This way we account for the screening of the nucleus by 
the electrons of the ion (note that this expression is still 
exact). Now we can take a limit T\ 3> 1 and apply mul- 
tipolar expansion to evaluate the integral in Eq. (10b I. 



The zeroth multipole cancels the first term in the right 
hand side of Eq. (Ill and we get 



h(n) 



N - 1 



- e. 



fi(n) 



7*(r0 p" k -\r,r)rr 2 dr. (12) 



Here we again omitted the angular factors and assumed 
that the sum runs only over states k which satisfy dipole 
selection rule. 

Let us assume that the solution of the homogeneous 
equation is localized at the distances r$. Then, for the 
larger distances n 3> T"j we can neglect the solution of 
the homogeneous equation lj6]) and write: 



/i( r l)ki>r 
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This expression can be simplified further if we note that 
at large distances the resolvent is approaching a constant 

N-l 
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In Eqs. (12 - 14) we did not make any assumptions 



about the wavefunction of the ion. In the single deter- 
minant approximation the integral in ( 14 ) is reduced to 



the sum over occupied orbitals of the opposite parity, 
PiPk = -1 0: 
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and we return to the Hartree-Fock case. 



operator: 



> 



Both equations (j 1 4|) and (|15| clearly lead to the asymp- 
tote (Q. For r , 
sum and 



the outermost orbital dominates the 



fi(r)\ r 



- v, fv(r). 



(16) 



Note that the power v\ can be larger than 2 if there is no 
dipole matrix element {fv\r\fi). In this case higher terms 
of the multipolar expansion are required. Expression Q 
shows that instead of the one electron binding energy 
of the valence orbital —s v one should use the ionization 
potential of the atom I N : 



r i exp 



(-V2W', 



(17) 



We see that Eq. (17) has much wider applicability than 



the Hartree-Fock approximation. On the other hand it 
is valid only for the distances r 3> r v , where for the 
neutral atom r v ~ 1. In fact, for the inner orbitals the 
exchange interaction starts to dominate over the direct 
Coulomb interaction much earlier, at r 3> r%, where r% is 
the radius of the inner orbital in question [3l [TTj . The 
more general form of Eq. (17) is given by Eq. ( |14| . In 
particular we can use them to estimate inner orbitals in 
the classically forbidden region. At large distances in the 
many electron atom all inner orbitals are asymptotically 
proportional to the outermost orbital up to a power of 
the radius. For shorter distances different terms of the 



sum (15) will dominate. 



Note that the system ( 10 ) is particularly useful for the 
highly charged ions, where 4 expansion is applicable. 
In this case the right hand side in Eq. (10a) is of the 



order and we can solve this system iteratively. For the 
first order corrections to the amplitudes fi we can use 
zero order wavefunctions ^ to calculate functions 

Wi^k- This way we can find analytical form of the first 
order corrections in 4 to the amplitudes fi. One simple 
example is considered in Appendix [B| 

In the end of this section we can say that asymptotic 
behaviour of the inner orbitals is changed by the entan- 
glement induced by the (anti) symmetrization postulate. 
This effect appears in the Hartree-Fock approximation, 
but survives for the correlated many-electron atoms. 



Observability of the one-particle asymptotics of the 
bound many-electron wave function 

Here we try to formalize the notion of the one-particle 
asymptotics of the bound many-electron wave function. 
In quantum mechanics one need to associate an opera- 
tor with any observable. When we discuss asymptotic 
behavior of the bound orbitals we mean that electron is 
registered at a given distance R from the origin and we 
simultaneously register the ion in a state with the hole 
in a given shell. This can be described by the following 
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and so on for the operators T 2 ^, T 3 ^, etc. When this 
operator is applied to the antisymmetric wave function 
fy N , we get: 



Expectation value of the operator Tf 1 for state 



(19) 



N 
4ttR 2 



(Mr)\S(r-R)\Mr)) . (20) 



Thus, we can say that functions fi indeed play the role of 
the orbital, whose long-range behavior we want to study. 

The observable ( 18 ) does not commute with the Hamil- 
tonian. The measurement of T t requires significant en- 
ergy. In order to detect an electron at a particular dis- 
tance R we introduce uncertainty in its momentum. That 
means that we interact with the whole many-electron sys- 
tem. Because of that the energy of the system is not con- 
served and we detect final ion with the energies E^ -1 , 
which are larger than the initial energy of the system 
E . The energy we need to detect position of the bound 
electron is E ^> £j. 

We see that the measurement of the observable Tj 
changes the energy of the system. This is different from 
the observations in the scattering theory where we detect 
particles at infinity. Such measurements can be done with 
arbitrary small momentum and energy transfer. Because 
of that in the scattering theory the energy of the system 
is conserved. These examples suggest that when during 
the measurement the energy of the system is conserved 
the exponent of the asymptote is given by the energies 
of the initial and the final state. On the contrary, when 
the measurement requires energy, the asymptote can be 
described by many exponents. 



Positron annihilation on inner electrons 



It seems that asymptotics Q and (17) are so differ- 
ent that it should be easy to prove experimentally which 
of them is correct. However these expressions coincide 
for the outermost electron shell, which gives dominant 
contribution to any physical processes at large distances. 
The inner shell asymptotic amplitudes are strongly sup- 
pressed. Though predictions for the inner shell contribu- 
tion to such processes from Eqs. ([2| and (17) can differ 
by many orders of magnitude, they still can be too small 
to be experimentally observable. 

Recently Amusia p~2] argued that asymptotics (17) 



should lead to an observable field ionization (FI) from 
the inner shells in the strong electric field F = Fz. He 
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estimated that the inner shell contribution scales as F 4 
and for the field F ~ 1 a.u. = 5 x 10 9 V/cm it can be 
on the order of 10~ 5 . Unfortunately such static fields are 
not achievable. If we use the low frequency laser field 
instead of the dc field, it will be hardly possible to dis- 
entangle tunneling FI from the multi photon processes. 
Therefore, it is unlikely that FI from the inner shells can 



be detected even for the asymptotics given by Eq. (17 1 



Here we want to draw attention to the positron anni- 
hilation on the inner electron shells [TBHTB] as a potential 
test of the expression ( 16 1. In the non-relativistic approx- 



imation the annihilation vertex is proportional to the 5- 
function. Thus, if we can control the final state of the 
ion, the annihilation cross section can be approximately 
linked to the observable T(R) defined by Eq. (18 1. 

In the typical experiment only one of the two annihi- 
lation gamma quanta is detected (see, e.g. the review 
[17|). For a low energy positrons an observed linewidth 
is determined by a Doppler broadening associated with 
the average momentum of the bound electron. Because 
of that the annihilation linewidth is given by |15j : 



I'M VN™eC 2 W3.7X 



|KeV. 



(21) 



Thus, the annihilation on the inner shells contributes to 
the wings of the line, while annihilation on the outermost 
shell gives the central peak. Consequently the accurate 
study of the annihilation line shape provides information 
about the inner shell contribution [13]. The inner shell 
annihilation was also detected directly with the coinci- 
dence technique for 7 quanta and Auger electrons 

Let us discuss Eq. (21) in more detail. The line shape 



of the annihilation 7 line is determined by the Fourier 
transform of the atomic orbital. It is clear that inner 
orbitals have wider spectrum and contribute to the wings 
of the line. However it is not so clear that the asymptotic 
part of the inner orbitals also contributes to the wings. 
Let us write the inner orbital as 



+ ■ 



(22) 



where </>*■ is a solution in a local potential and 4>\ describes 
exchange induced asymptotic tail. Let us write <pf as 



C 



;4>v(r) ■ 



(23) 



This function has correct asymptotic behaviour ( 16 ) with 
Vi = 2 and parameter a is introduced to insure that at 



short distances 



> 



The Fourier transform of this 



function is the convolution of the transforms of the mul- 
tipliers. The width of the convolution is the sum of the 



widths of the components, i.e. Pf 



P„, where 



P,, 



is the width of the Fourier transform of the va- 



lence orbital <p v . The exchange term starts to dominate 
the tail at the distances between r* and r v [2] . Thus, the 
cutoff parameter a should be chosen from the interval 
r*j < a < r„. Even if we take a close to r v we get the 
width two times larger than for the outermost orbital. 



However, a is closer to r, since it plays a role of the cut- 
off factor in the expansion over r < /r > in the exchange 
Coulomb interaction between inner and outer electron. 
For a ~ r, we return to Eq. (21 1. In both cases the an- 



nihilation on the asymptotic part of the inner orbital (f>f 
indeed contribute to the wing of the 7 line. 

We conclude that in the annihilation experiment we 
can control the final state of the ion. Unfortunately it is 
not the case for the position of the annihilated electron. 
Still, the positron can not penetrate deep into the atom 
because of the strong repulsion from the nucleus. There- 
fore, the annihilation should predominantly take place at 
the edge of the atom. 

Let us estimate classical turning point for the thermal 
positron with the energy e p ~ 300 K ~ 10~ 3 a.u. Atomic 
potential seen by the slow positron has the form [18 : 



U(r) = Q{r)/r-a/2r 4 



(24) 



where a is static polarizability of the atom. At large 
distances we can parameterize effective charge Q(r) as 



N, 



Q(r) = exp (-2V=2^(r - r v )) , 



(25) 



where e v and N v are the electron energy and the number 
of electrons for the outermost atomic shell. The radius 
r v can be defined so that Q(r v ) = -g-, i.e. the screening 
of the nucleus at the distance r v by the outer shell is 
reduced by 50%. For the neutral closed shell atom r v ~ 1 
and N v = Al v + 2. The classical turning point for the 
positron r t is given by the equation: U(rt) = £ p <C 1, or 
neglecting positron energy e p : 



1 N r A 

r t~r v + In 

2|2e t ,| i /^ a 



r " + 2|2^p' (26) 



where we approximated the logarithm by unity. We see 
that rt is only about an atomic unit larger than the outer 
shell radius r v . This is much smaller than the typical tun- 
neling distance in the dc electric field, making positron 
annihilation much more sensitive to the inner shell con- 
tributions. On the other hand, the distance r t is much 
larger than the inner shell radii, so annihilation on the in- 
ner electrons should depend on the asymptotic behaviour 
of the inner orbitals in the classically forbidden region. 
The estimate ( 26 ) for the positron turning point r t is 



close to r v where Eq. (17) is not applicable, but Eq. (16) 
should hold. 

As a first example we considered Be atom in different 
configurations and compared asymptotics of the Is or- 
bital. We used relativistic Hartree-Fock-Dirac code [H] 
for three closed shell configurations: Is 2 , ls 2 2s 2 , and 
ls 2 2p 2 y 2 . Fig. jlj shows the ratio of the Is orbitals for all 

three cases. One can see that for the configurations Is 2 
and ls 2 2s 2 the long distance behaviour is similar. The 
small difference is caused by a 20% change in the Is en- 
ergy. We conclude that here numerical results agree with 
Eq. (I2J). The Is energy for the configuration ls 2 2p 2 ^ 9 lies 
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FIG. 1: Ratios of the Is and 2p orbitals for different configu- 
rations of Be. Vertical arrows show rms radii of Is and 2pi/ 2 
orbitals and the positron turning point rt- 



FIG. 2: Ratios of the 5p 3/2 , 5si/ 2 , 4d 5/2 , 4p 3/2 , and 4si/ 2 
orbitals in Xe atom. Vertical arrows show rms radius of 5p 3 / 2 
orbital and the positron turning point r t . 



between the values for two other configurations, but the 
asymptotics is absolutely different in agreement with Eq. 
Q. This is in consent with the statement that Eq. ([!]) 
applies for the systems with occupied shells with I ^ [1] . 
Fig. [T] also shows the ratio of the 2p 1 / 2 and Is orbitals 
for the configuration ls 2 2py 2 . We see that this ratio 

first grows exponentially, but then stabilizes at r > 3 in 
agreement with (16). 



Figure [T] shows that exchange interaction start to de- 
termine behaviour of the Is wave function at the dis- 
tances r < r 2p , i.e. near the main maximum of the 2p 
orbital. At such intermediate distances, which lie far be- 
hind the classical turning point for the inner electron but 
in the localization domain of the outer electrons, the ex- 
change interaction is suppressed only by a power of the 
radius HE]. 

For the Be atom in ls 2 2p^ 2 configuration the estimate 



(261 gives r t s» 4.2. Thus, the distances where annihi- 
lation can take place are r ~ 4. At such distances the 
difference between Is orbitals with and without exchange 
interaction (configurations Is 2 and Is 2p\, 2 respectively) 
is about one order of magnitude. Therefore, the exchange 
assisted tunneling can enhance annihilation rate (which 
is proportional to the probability density) by approxi- 
mately two orders of magnitude. As we see from Fig. [T] 
the ratio of the 2p and Is orbitals at these distances is on 
the order of 10 2 . Therefore, for the Be atom in the con- 
figuration ls 2 2p 2 we can expect the Is contribution to 
the annihilation rate to be on the order of 10 -4 , instead 
of 10 -6 without exchange assisted tunneling. 



Let us consider now a more realistic example. The 
inner shell contributions to the annihilation 7 line was 
studied by Iwata et al. [T3] for several noble gases. For 
Xe the total probability of the annihilation on the inner 
shells id, 4p, and 4s was found to be 2.4%. On Fig. [2] we 
plot the ratios of the orbitals for Xe at the distances from 
1 to 8 Bohr radii, where these orbitals do not oscillate. 



For Xe the estimate ([26) for the classical turning point 
gives r t ~ 3.1. At such distances the ratio of the 5p and 
4c? orbitals is close to 10. Taking into account that d shell 
has 10 electrons while p shell has only 6, we conclude that 
2% contribution of the 4c? shell agrees with our simple 
consideration. The rms radius for the 4d orbital in Xe 
is 0.95 a.u., so such a large contribution is due to the 
asymptotics ( 16 ). 



Correlations play very important role in the process 
of the positron annihilation. They lead to the huge en- 
hancement of the cross section. The electron core polar- 
ization leads to the attractive potential for the positron 
(24) increasing the annihilation rate. Usually this ef- 



fect has the size of a typical correlation correction for 
many electron atoms. The dominant correlation correc- 
tions come from the positron-electron correlations [18) . 
This type of correlations can lead to the virtual forma- 
tion of the positronium and neutralization of the positron 
charge. As a result, the positron can penetrate deeper 
inside the atom. As we saw above, the classical turn- 
ing point for the positron r t lies in the region where 
the electron density of the atom is exponentially decreas- 
ing, so even small correlation corrections to r t can dras- 
tically increase the annihilation rate. However, it was 
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shown, by Green and Gribakin [20 that these correla- 
tions equally affect outer and inner shell contributions 
and only slightly change the annihilation line shape. Note 
that this result indirectly confirms that the densities of 
the inner and outer electron in the annihilation region 
are almost proportional to each other, as suggested by 
Eq. JTl. 



In this section we considered direct annihilation on the 
inner shell electron. We argued that this process is en- 
hanced by the exchange interaction between inner and 
outer shells. The inner shell hole can be also formed in 
the annihilation on the outer shell followed by ionization 
of the inner shell caused by the changed atomic potential. 
In the lowest order of MBPT this process is described by 
the Coulomb interaction between the inner and the outer 
shell holes. This Coulomb integral is the same as the ex- 
change integral considered above with substitution of the 
one final orbital by the orbital in the continuum. Such 
processes were considered by Dunlop and Gribakin [15] . 



Conclusions 

In this paper we argue that the exchange assisted tun- 
neling of the electrons from the inner atomic shells is not 
an artifact of the Hartree-Fock approximation, but an 
observable effect, at least for the intermediate distances 
r > 1. In particular, it can be observed in the annihila- 
tion of the slow positrons on the many electron atoms. 
The annihilation on the inner shell electrons forms the 
shoulders of the experimentally observed 7 lines for the 
noble gases [13]. It can be also directly detected in the 
coincidence experiments where 7 quanta are registered si- 
multaneously with the Auger electrons [T3] . This process 
takes place at the distances, comparable to the size of the 
outermost atomic shell. Because of that it is much less 
suppressed than the inner shell dc field ionization, which 
may be difficult to observe. On the other hand, without 
the exchange assisted tunneling of the inner electrons to 
this region, the inner shell annihilation would be unob- 
servably small. At present we do not see any realistic 
experiment to test the inner electron asymptotics at the 
large distances r > 1 in atomic physics. However such 
asymptotics can be important in the condensed matter 
physics [3J. 
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FIG. 3: Double well potential with two localized and one 
delocalized states. 



Appendix A: Double well potential 

Let us consider a double well potential. We assume 
that there are two localized states L and R and one 
delocalized state D with energies sl, £r, and Ed (see 
Fig. [3]). If El = £r = £lr, then the true eigenstates 
are ip± = (i/'l±'0-r)/'s/2 and there is exponentially small 
energy splitting between them: 



£± =£lrT$, 



(Al) 



The one-particle spectrum of this system is e + , e_, and 
£d- 

If we have two non-interacting electrons in this poten- 
tial in a triplet state, the spectrum still includes three 
levels: 



E„ 



2e 



LR ■ 



£lr + £d T S ■ 



(A2) 



If we switch on interaction l/\x\ — x 2 \ between the elec- 
trons, we get additional splitting between levels Ef, and 
E c from the exchange integral: 



E c -E b = 2(5 + A) , 

A = ^L(x 1 )ip D (x 2 ) 
2 



\X!-X 2 \ 



(4>l\x\iPd}{iPd\x\'<Pr) 



(A3) 



(A4) 



J LR 



Here we left the first non-zero term of the multipolar 
expansion (xlr. is the distance between the minima of 
the potential). 

We see that for the non-interacting electrons the split- 
ting is caused by the tunneling through the barrier and is 
of the order of the overlap integral (V'lIV'-r)? which expo- 
nentially depends on the distance xlr. The exchange in- 
teraction with the delocalized electron results in the split- 
ting which depends on t his distance as x~£ R . The dipole 



matrix elements in (A4) are not suppressed because the 



state D has large overlaps with both localizes states. We 
conclude that the exchange interaction can lead to the 
long range interaction between localized states. 
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Appendix B: Z 1 expansion for the He-like ion 



The system ( 10 ) can be used to form the i expansion 



of the eigenfunction of the highly charged ion. Let us 
consider He-like ion with Z 3> 1. Making substitution 
£j = Zri we write Hamiltonian ^ as: 



Z 2 Hi 2) 



= Z Z 



1. 1 1. 1 



1 1 



In the new variables we have 

g (2)^(2) ; e (2) _ ^-2^(2) 



^(2)^(2) 



(Bl) 



(B2) 



We present the solution in the form 
system (10 1 where 



and need to solve 



1 



.(2) 



(B3) 



Here the functions </> ni ,Zi are hydrogenic. Explicit small- 
ness in ( j B3 h and, therefore, in the right hand side of 
Eq. ( lOab allows expansion in Z^ 1 . 



Let us consider the triplet state Is2p m —Q. The zero- 
order orbital functions have the form: 



*o 2) = -^(^(1)^.(2) - <M1)<W2)) , (B4) 



1 



fopo (0 



e 



V47T 

cos 9 
Ay/21: 



-f/2 



(B5) 
(B6) 



Note that already in the first order in Z^ 1 we have in- 
finite number of channels. However, these channels cor- 
respond to the different final states of the ion and can be 
considered independently. We are interested only in the 
first order corrections to the orbitals f\ and f 2 , which 
correspond to the ion either in the Is, or in the 2p state: 



fi - 4>2 Pa + fill > 

h = 4>u + /£> + /i 1} 



(B7) 
(B8) 



Note that angular and parity selection rules allow mixing 
of s- and d- waves for the amplitude f 2 - However, we will 
be interested only in the correction to the s-wave . 



For the orbitals (|B5| and (B6) the functions Wik can 



be found analytically using Eq. ( B3 ) : 



Wi,i 
W 1>2 

W 2 ,i 

W 2 . 2 



2cos0 
~243Z£ 2 

+ 72£ 2 + 96£ + 64) 

2 



(64-e- 3?/2 (27C 3 



( 6 4-e~ 3 «/ 2 (27£ 3 



729Z£ 2 cos 
+ 72£ 2 + 96£ + 64) 



V2Z£\ 24 6 « 



6£ 2 + 18^ + 24) 



(B9) 



(BIO) 



(Bll) 



(B12) 



In Eq. (Bill we left only the term for which the prod- 



uct 4>2p W2,i has no angular dependence and, th eref ore, 
contributes to the s-wave part of the amplitude ( B8 ) . 



The first order corrections and satisfy equa- 
tions: 



1 1 



~2 A ' £ + 8 ) 



»2p 



Wi 



b ls W h2 , (B13) 



IJ 2po 



w 2 



(B14) 



These equations are valid at all distances. We see that 
the right hand sides of both of them have three differ- 
ent exponents, — £/2, — £, and — 5£/2. The long-range 
asymptotics depends on the weakest exponent — £/2, 
which corresponds to the binding energy of the electron 
2p. We conclude that the first term of the Z~ x expan- 
sion has asymptotics in agreement with the general case 
discussed above. 
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